Thermal robustness of multipartite entanglement of the 1-D spin 1/2 XY model 
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We study the robustness of multipartite entanglement of the ground state of the one-dimensional 
spin 1/2 XY model with a transverse magnetic field in the presence of thermal excitations, by 
investigating a threshold temperature, below which the thermal state is guaranteed to be entan- 
gled. We obtain the threshold temperature based on the geometric measure of entanglement of the 
ground state. The threshold temperature reflects three characteristic lines in the phase diagram of 
the correlation function. Our approach reveals a region where multipartite entanglement at zero 
temperature is high but is thermally fragile, and another region where multipartite entanglement 
at zero temperature is low but is thermally robust. 



I. INTRODUCTION 

In quantum information, we manipulate quantum sys- 
tems using unitary operations and measurements. To 
achieve such manipulations with high accuracy, well- 
controllable quantum systems are required. Physically, 
well-controllable quantum systems are likely to be de- 
scribed by microscopically controllable Hamiltonians. 
We may describe such systems as artificial systems. On 
the other hand, in natural physical systems, it is rare 
to have a microscopically fine-controllable Hamiltonian. 
Typically we cannot switch on and off individual interac- 
tions on demand and it is more likely that we only have 
control of a few macroscopic parameters of the Hamilto- 
nian, such as the a magnetic field. 

Entanglement is a non-local correlation that does not 
exist in classical mechanics, and it is considered to be 
a key for understanding the power of quantum compu- 
tation and quantum protocols [T], when comparing to 
their classical counterparts. Although the controllability 
of systems is different for artificial systems and for natu- 
ral physical systems, entanglement is an indicator of the 
quantum signature for both systems. Recently, investi- 
gations of quantum effects from the view point of entan- 
glement in natural physical systems attracts increasing 
interest 0. 

One interesting result on the relationship between en- 
tanglement and quantum properties of natural physical 
systems regards entanglement of their ground state and 
quantum phase transitions (QPTs) [3J. QPTs are phase 
transitions at zero temperature induced by varying pa- 
rameters in a Hamiltonian [4 . It is numerically shown 
in many spin models that measures of bipartite entangle- 
ment of the ground state, for example, concurrence [5] of 



the reduced two-spin state of the ground state, or block 
entropy, has non-analyticity in the vicinity of the QPT 
points [3J. 

Since isolating the systems completely from the envi- 
ronment is not easy for natural physical systems where 
we do not have fine control of microscopic variables, it is 
difficult to achieve zero temperature. Thus, it is interest- 
ing to understand the relationship between entanglement 
properties and quantum effects corresponding to QPT of 
natural physical systems at non-zero temperature. There 
have been several works investigating bipartite entangle- 
ment of non-zero thermal states of natural physical sys- 
tems [B], and for many models, it is found that concur- 
rence and logarithmic negativity of thermal states are 
nonzero at finite temperature. 

However, there is no guarantee that such bipartite en- 
tanglement fully describes the properties of the QPT. In 
general, eigenstates are fully entangled in a multipartite 
manner. From the study of multipartite entanglement 
in quantum information, it is known that some entangle- 
ment properties (e.g., SLOCC convertibility [7]) of mul- 
tipartite states cannot be characterized only by bipar- 
tite entanglement even for pure states. For mixed states, 
there even exists a state where there is no bipartite entan- 
glement between any two division, but still contains mul- 
tipartite entanglement [8]. It has also been shown that 
resonating-valence-bond (RVB) states, which are used to 
describe many natural many-body systems [3], contain 
multipartite entanglement but have a negligible amount 
of two-site entanglement [10] . 

This has led to investigations using multipartite en- 
tanglement, though this avenue is much less well devel- 
oped, and there are fewer results. The main obstacle 
is the difficulty in calculating many of the multipartite 
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entanglement measures for large physical systems. Due 
to this, one approach has been to consider multipartite 
entanglement witnesses, which are generally much more 
mathematically tractable (see e.g. jTTJ Q2]). This has 
given bounds on the minimum energy allowed for ther- 
mal states of the spin 1/2 Ising model and XX Heisenberg 
in a transverse magnetic field for example |T3} [H] ■ Also 
in this direction, in |15j . the authors constructed a multi- 
partite entanglement witness for mixed systems including 
thermal stabilizer states related to (multiple-point) cor- 
relation functions. Another direction has been to look 
for entanglement measures which arc easier to calculate. 
In [16 j the geometric measure of entanglement is studied 
for the ground state for the ID spin 1/2 XY model over 
QPTs, and it is seen that over a phase transition, the 
geometric measure behaves singularly. But in general, 
due to the computational difficulty of analyzing mixed 
multipartite entanglement, multipartite entanglement in 
the thermal states of natural physical systems and its 
relationship to QPT-like behavior of the system has not 
been well understood. 



In this paper, we study multipartite entanglement of 
a thermal state in the ID spin 1/2 XY model with a 
transverse magnetic field to investigate the relationship 
between multipartite entanglement and the QPT-like be- 
havior of natural physical systems at non-zero temper- 
ature. Instead of analyzing multipartite entanglement 
of an exact thermal state of the system, we investigate 
how thermal excitations affect the entanglement of the 
ground state. We employ the threshold temperature, 
which is the temperature below which a thermal state 
is certainly entangled, defined in Ref. 17J. We derive the 
threshold temperature of multipartite entanglement of 
the ground state measured by using a geometric measure 
[18j and information of population of the ground state. 
This threshold temperature is an indicator to present how 
robust multipartite entanglement of the ground state is 
in the presence of thermal excitations. In addition, we see 
that this approach extends the identification of borders 
by entanglement seen in [TB] , showing that the threshold 
temperature varies singularly at all of the borders of the 
phase diagram of the system. 



This paper is composed as follows. In Section [IT] we 
present a threshold temperature and show how an energy 
gap between the ground energy and the first excited en- 
ergy affects to the threshold temperature. We describe 
the ID spin 1/2 XY model with a transverse magnetic 
field and show a phase diagram and corresponding multi- 
partite entanglement at zero-temperature in Section |III| 
We present our results and analysis of the threshold tem- 



II. THRESHOLD TEMPERATURE 

A. Definition 

We consider systems described by a Hamiltonian H 
and denote the eigenenergies and eigenstates as {Ei} 
and {|e7)} 7 =i,--- ,di , respectively, where i — 0,1,- •• and 
di is the degree of degeneracy of Ei. By representing 
each eigen-subspace for Ei by using a projector = 



E 



di 
7=1 



(e?\, each degenerate mixed state correspond- 



ing to Ei is written by pi — Pi/di. If the system is 
in the thermal equilibrium at temperature T, the ther- 
mal state is written by p t h{T) = exp[-H/(k B T)]/Z(T) 
where k B is the Boltzmann constant and Z(T) — 
J2i exp[— Ei/(ksT)] is a partition function of the system. 

For multipartite systems, calculating eigenenergies and 
eigenstates is computationally difficult in general. Fur- 
ther, even for the cases that eigenenergies and eigenstates 
are somehow obtained and we can explicitly write down 
the thermal state pth(T), finding out whether the given 
thermal state pth(T) is separable or not is another com- 
putationally very difficult problem, since no simple sepa- 
rability criterion is known for multipartite mixed states. 

Instead of using all the information of eigenenergies 
and eigen-subspaces of the Hamiltonian, we analyze this 
system by just taking into account the information of 
the ground state po and its population d^po, where 
p (T) = exp [-E /(k B T)]/Z(T) is the population of the 
ground state, to investigate multipartite entanglement of 
a state under thermal excitations. Using this limited in- 
formation, we will find bounds on the minimum temper- 
ature for which entanglement can be guaranteed. This 
is the concept of a threshold temperature introduced by 

El. 

To define the threshold temperature, we first use a 
distance-like measure of multipartite entanglement, the 
robustness of entanglement [19|. The robustness of en- 
tanglement R(p) for a general state p is defined by 



R(p) := min{i 



1 

TTt 



(p + tuj) G Sep} 



(1) 



where Sep denotes a set of all separable (mixed) states. 
If we consider (p + tu>)/(l + t) as a mixture of a state p 
with another state ui with the their populations 1/(1 + t) 
and t/(l + t), respectively, we can interpret 1/(1 + R(p)) 
as the population of the state p when lo is a state which 
destroys entanglement of p most effectively. 

The thermal state p t h (T) can be interpreted as a mix- 
ture of the ground state po and the-rest state ui' as 



Pth{T) = d pa(T)p + (1 - d Q p (T))u)' 



(2) 



where w' = 1/(1 - d p (T)) J2i=iPi( T )P*- % definition 
([I]), the thermal state pth(T) should be entangled at the 
temperature T if the corresponding populations satisfy 
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Po(T) > 



d (l + R( Po ))' 



(3) 
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Since the state u>' is not necessarily the state u> that de- 
stroys entanglement of po most effectively, this is only a 
sufficient, not necessary condition for entanglement. 

In general, evaluation of the robustness of entangle- 
ment is still a computationally difficult problem since it 
requires searching the mixed state lo over all the Hilbcrt 
space. However, it has been shown [TS] that a lower 
bound to the robustness of entanglement is given by the 
geometric measure G(p) [5111 |2"T] 

d[l + R(^)]>2 G ^\ (4) 
a 

where P is the support of the state p and d is the rank 
of the support P. The geometric measure is defined by 

G(p) :=-log 2 [ max (3|p|3>] (5) 

|*)GPro 

where Pro denotes a set of all product pure states. The 
search of the pure product state |$) greatly reduces the 
required computational efforts comparing the search of 
the mixed state u> for R(p). 

Although the geometric measure is a measure of en- 
tanglement for pure states, it is not for mixed states. 
However, it still gives a nontrivial lower bound of the 
robustness of entanglement |21j . Therefore, we define a 
lower bound threshold temperature T t u using the geo- 
metric measure which satisfies 

Po (T th ) = 2- G{ ^\ (6) 

We stress that this is a temperature below which the ther- 
mal state is certainly entangled. It is a lower bound of a 
temperature for the state to be entangled and a thermal 
state above the threshold temperature is not necessar- 
ily separable. In this sense it fills the same role as an 
entanglement witness. 

B. Effects of the first energy gap 

Before investigating the threshold temperature T t h of 
the 1-D spin 1/2 XY model, we analyze how an energy 
gap between the ground energy and the first excited en- 
ergy affects to the threshold temperature. It is intuitively 
expected that a system with a smaller energy gap has a 
lower threshold temperature, since the population of the 
ground state decreases quickly by increasing temperature 
in such cases. 

To clarify this intuition,we take into account the popu- 
lation of the first excited state. Our resulting relationship 
is that 

(l- e -4r)p (T) > - rp * , =*> Pth(T) is entangled, 

; (7) 
where A = E\ — Eq. The proof is given in 
Appendix [A] By using the inequality Q, we de- 
fine a gapped threshold temperature T t ^(A) by (1 — 



This gapped threshold temperature T t ^(A) is always 
smaller than the corresponding threshold temperature 
T t h. The reason of T t ^(A) < T t h is that a bound of 
the robustness of entanglement to obtain the relation- 
ship ([7]) is too loose (see Appendix |A| . We need to 
be aware that this does not imply that the threshold 
temperature decreases if we take into account the first 
energy gap. The gapped threshold temperature is use- 
ful to understand how the gapped threshold temperature 
is affected by varying the first energy gap although the 
gapped threshold temperature is not much use to evalu- 
ate a bound of the threshold temperature, given by Eq. 
([5]). In the rest of this section, we examine the gapped 
threshold temperature to study the effect of the first en- 
ergy gap. 

Since (1 - e- A /( feBT ))po(7 1 ) > (1 - e - A '/( fe ^)) Po (T) 
for any T and any A > A', and the geometric measure 
of the ground state is independent of the energy gap A, 
we can conclude that T t h(A) > Ttfe(A') for finite gaps 
A > A'. Thus, a smaller energy gap results in the lower 
gapped threshold temperature as we have expected. We 
present an example to display this relationship. 

We consider a m-partite spin 1/2 system. Suppose that 
the ground state is a generalized W-state given by 



I^H^aU-'-D + IIT •••4} + ---+IU---T» (8) 



and the population of the ground state is given by 
Po(T) = ^"Ifl^^T) where S is a constant and D = 2 m 
is the whole dimension of the system. We assume that the 
first energy gap is kS. Since it has been shown that the 
equality in the inequality Q holds as an equality for the 
W-state by [18] and its geometric measure is calculated 
to be (m - l)log 2 ^ gD] [21], the gapped threshold 

temperature T t ^(«;5), is given by the root of the equa- 
tion 

(1 _ e -«<5/(fc B T th )^ 1 _ e -8/(k B T th )^ = ^ ^ 

in the thermodynamic limit (m — > 00). 

We vary the first energy gap kS to analyze the gapped 
threshold temperature. Note that we fix the population 
of the ground state by adjusting other eigenenergies. We 
present the gapped threshold temperature T t h(nS) as a 
function of the k in Fig. [I] We can see that the gapped 
threshold temperature changes drastically with n in the 
region of small k although it does not change much for 
large k. This example shows that the existence of the first 
energy gap A, namely whether A = or A ^ 0, is very 
influential on the gapped threshold temperature since the 
gapped threshold temperature changes drastically with A 
if the first energy gap is small. 
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FIG. 1: (Color online) The threshold temperature Tthi^S) as 
a function of the k. We set Boltzmann constant fcg as one 
for simplicity. The dashed line is the threshold temperature 
T t h defined by Eq. (J5|. By definition, the gapped threshold 
temperature T t h(n.8) converges to the threshold temperature 
Tth in the limit of k — > oo 



III. ID SPIN 1/2 XY MODEL WITH A 
TRANSVERSE MAGNETIC FIELD 

The dimensionless Hamiltonian of the ID spin 1/2 XY 
model with a transverse magnetic field is given by |22] 

H XY = - 2^y-^ a 3 a J+i + ~2~ a J a J+i + ha J J' 

3=1 

(10) 

where er" (a — X, Y, Z) are Pauli matrices acting on a 
spin at j-th site, r\ £ [—1, 1] is a parameter of anisotropy 
and h > is an external magnetic field. We consider 
that 77 and h are normalized by the coupling energy of 
the interaction term in the XY model Hamiltonian. For 
convenience, we consider only a system with even number 
of particles N here, but the following arguments can be 
similarly applied to the case of odd N. We also employ a 
periodic boundary condition cr^ +1 := erf (a — X,Y,Z). 

Although this XY model has been much studied in 
condensed matter physics, it is not so familiar in quantum 
information. We present the phase diagram of the XY 
model at zero temperature and the entanglement of the 
ground state. A brief explanation of the eigenstates and 
the eigenvalues is in Appendix |B) 



A. Phase diagram at zero temperature 

We present the phase diagram of the XY model char- 
acterized by two point correlation functions at zero tem- 
perature as in Fig (2] In different phases, two point 
correlation functions C a (r) between two points sepa- 
rated by distance r, defined by C a (r) := (ct"<t" +1 ) for 
a = X,Y, Z , have different characteristics. Although this 
two point correlation functions generally decrease expo- 
nentially with the distance r, C a (r) ~ exp[— |] where £ 
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FIG. 2: (Color online) The phase diagram at zero temperature 
in the XY model. These phases are characterized by two 
point correlation functions. The ground states degenerate 
in the shaded region. The dashed line (rj = and h > 1) 
in the phase diagram is not a phase transition line but the 
correlation function on this line is constant with respect to 
the distance r similar to the C-line. 

is a constant called a correlation length, the rate of de- 
crease is slower and is generally polynomial with r on a 
border of phases. 

In the phase diagram of the XY model at zero tem- 
perature, there are three phases, phase 1, phase 2 and 
phase 3, which are divided by three borders referred to 
as a vertical line (V-line), a horizontal line (H-line) and 
a constant correlation line (C-line). The phase 1 is the 
region of < h < \J\ — rj 2 , the phase 2 is the region of 
\J\ — T] 2 < h < 1 and the phase 3 is the region of 1 < h 
(see Fig. [2}. 

The correlation functions have been calculated for 
large r in the thermodynamic limit |23j . In the phase 
1, the correlation functions oscillate with the distance 
r with exponentially decreasing envelope. In the other 
regions, there is no oscillation. The rate of decrease in 
the phase 2 is different from that in the phase 3. Then 
they are considered as three different phases. We describe 
properties of the system on the V-line (h — 1), the H-linc 
(77 = and < h < 1) and on the C-line (h 2 + if = 1). 



1. V-line and H-line 

On the V-line and on the H-line, a quantum phase 
transition (QPT) occurs. QPTs are defined by the non- 
analyticity of the ground energy with respect to param- 
eters in a Hamiltonian, and a line on which a n-th order 
QPT occurs is called a n-th order QPT line. The V-line 
is a second order QPT line since the second derivative 
of the ground energy diverges. Due to this second or- 



der QPT, the correlation functions C a (r) decrease only 
polynomially with the distance r |23j . 

The H-line is a first order QPT line. The first deriva- 
tive of the ground state energy diverges and there is a 
crossover of the eigenstates on the line for the ground 
state. This first order QPT also results in the polyno- 
mial decreasing of the correlation functions C a (r) with 
the distance r l24l. 



2. C-line 

The C-line is not a QPT line. However, the two point 
correlation functions C Q (r) are constant with respect to 
the distance r. In the thermodynamic limit, they are 
given by 



C x {r) 



2rj 
I + 77' 



C Y (r) = 0, 
C Z (r) 



(11) 

(12) 
(13) 
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FIG. 3: (Color online) The geometric measure of the ground 
state of the XY model (upper) and its 2D plot (lower). The 
geometric measure is numerically calculated for the system of 
N — 80. Note that the geometric measure in the region of 
r\ £ [—1,0) is the same as that in the region of r\ G (0, 1]. 
The same result had been originally obtained for the region 
of 7>G [0,1] inRef. [TBJ. 



where mz is the magnetization in the Z-direction defined 
by mz — (Eili a i /N) [23 . These constant correlation 
functions originate from the fact that the ground state on 
the C-line approaches to a separable state. The explicit 
form of the degenerated ground states on the C-line is 
given by |9+) = (cos f | f) + sin f | i))® N and |9_ 

(cos f I T) - sin f I I) )® N where cos 6 = 

Although these two ground states are both separable 
states, we note that the total degenerate ground state 
of this system not precisely separable. Since the ground 
states are on the C-line degenerate, the total ground state 
is a mixed state given by the projector onto the subspace 
spanned by two ground states |0±), which is written as 




1 



Po = ^(|e+>(e + | + |ei)(ei|), 



(14) 



where |Oi) is a state orthogonal to |0+). The state |9i) 
is given by 



|e±) 



1 



Vi - cos 2N e 



e_)-cos iV e|e+)). (15) 



The state |0:j:) is easily shown to be entangled. It is 
known |26j that a mixture of any pure entangled state 
with any separable pure state is entangled. Thus, the 
mixed ground state po is entangled. However, as we 
will see the next subsection, entanglement of the mixed 
ground state po is weak and this weak entanglement re- 
sults in the constant correlation functions. In addition, 
we emphasize that the ground energy itself has no char- 
acteristic on the C-line. The C-line is characterized only 
by properties of the ground state. 



B. Geometric measure of the ground state 

We present the geometric measure of the ground state 
of the XY model. The geometric measure had been orig- 
inally obtained in Rcf. [15] for the region of rj G [0,1], by 
considering the pure ground state (either |ej(/i, 77)) or 
|e^~(/\, 77)) defined in Appendix pb although the ground 
state for h 1 is a mixed state due to its degeneracy (see 
Appendix |B| . In our calculation, we consider the mixed 
ground state po for the region of h < 1 . It turns out that 
the geometric measure of the mixed ground state coin- 
cides with that of each pure degenerated ground states in 
the thermodynamic limit. This is because the scaling be- 
havior of the geometric measure of the two degenerated 
ground states are the same as in Ref. [16 . Our numerical 
calculation of the geometric measure of the ground state 
of the XY model is presented in Fig. [3] 

The first derivative of the geometric measure diverges 
on the V-line. This divergence is expected to originate 
from the fact that the ground state depends on the vari- 
ables, which have the non-analyticity on the V-line (see 
Appendix [B]). The non-analyticity of the variables results 
in the dramatic change of the ground state around the 
V-line (see Appendix [C] for more detail). On the H-line, 
the geometric measure is non-differentiable because of 
the same reason. We note that the system in the case of 
rj = has a different symmetry from other regions. 

The geometric measure on the C-line does not have 
distinguished characteristics to the same degree as on 
the V-line and the H-line. The geometric measure ac- 
tually takes its local minimum on the C-line with re- 
spect to parameters (h, 77), which results in constant cor- 
relation functions (see Subsection III A I. However, since 
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FIG. 4: (Color online) The threshold temperature in the XY 
model (upper) and its 2D plot (lower) . It is calculated in the 
system with N = 80. We set Boltzmann constant ks as one. 
We can clearly see the V-line, the H-line and the C-line. We 
note that the line of 7] = 0, h > 1 corresponds to the dashed 
line in the phase diagram (see Fig. [2}, which is not a QPT 
line. 



the geometric measure of the ground state is still small 
around the C-line and changes smoothly around the C- 
line. the C-line is difficult to identify (see the 2D plot in 
Fig. Ill 13 I . Note that the two point correlation functions 



are constant only on the C-line. 



IV. THRESHOLD TEMPERATURE IN THE XY 
MODEL 

A. Threshold temperature and the phase diagram 

Since the exact population of the ground state of the 
XY model is known J23], numerical calculations of the 
threshold temperature is straightforward by using the 
values of the geometric measure of the ground state. We 
present the result of our numerical calculations in Fig. 
g| First, we focus on the singularity of the threshold 
temperature and analyze the correspondence between the 
singularity of the threshold temperature and the phase 
diagram at zero temperature. 

On the V-line, the first derivative of the threshold tem- 
perature diverges. This is because the threshold temper- 
ature defined by po{h,r] : T t h) — 2~ G ( p °) inherits the 
divergence of the first derivative of the geometric mea- 
sure on the V-line and the first derivative of the pop- 
ulation of the ground state is finite. Similarly, on the 
H-linc, the threshold temperature is non-differentiable. 
Its non-differentiability originates from that of the geo- 
metric measure of the ground state on the line. 

The threshold temperature drops dramatically on the 
line of h > 1, 77 = and on the C-line. These character- 
istics originate from the local minimum of the geometric 







0.8 






0.6 






n a 










>S 


0.2 




1 





0.5 



FIG. 5: (Color online) The threshold temperature shown in 
Fig. [4] The (red) solid line with small squares, the (green) 
dashed line , the (blue) dotted line , and the (purple) solid line 
correspond to 77 = 0,0.4,0.9, 1.0, respectiv ely. Th e threshold 
temperature has the singularity at h = \J\ — r\ 2 (on the C- 
line) and at ft = 1 (on the H-line). 



measure on both lines. Since the ground state on the line 
of h > 1, rj = is separable, which also results in constant 
correlation functions, the geometric measure G(po) and 
the threshold temperature T t h is zero on the line. On the 
C-line, the geometric measure G(po) takes local minimum 
with respect to the parameters {h,rj), which leads to a 
local maximum of 2~ G ( p °K Since the population of the 
ground state po(h,r] : T) decreases monotonically with 
temperature, the local maximum of 2~ G ( p °) yields a kink 
of the robustness threshold temperature as long as the 
change of the population po(h, rj : T) by varying (h, rj) is 
less than that of the geometric measure G(po(h,r) : T)). 
Thus the threshold temperature can clearly identify the 
region where entanglement of the ground state is the 
weakest. The constant correlation functions are expected 
to originate from weak entanglement, the ground state in 
the region where the threshold temperature dramatically 
drops will give constant correlation functions. 

We note that the shape of the singularity of the thresh- 
old temperature depends on the characteristics of the 
line. On the first order QPT line, the threshold tem- 
perature takes a kink as we can see on the H-line. On 
the second order QPT line, the first derivative of the 
threshold temperature diverges, which means the thresh- 
old temperature dramatically changes by varying param- 
eters around the second order QPT line. We can see it 
on the V-linc. Finally, if the correlation functions are 
approximately constant on a line, the threshold temper- 
ature drops dramatically on the line as shown on the 
C-line and on the line of ft. > 1 and rj = 0. 



B. Threshold temperature in the thermodynamic 
limit 



We investigate the threshold temperature in the ther- 
modynamic limit in this subsection. In principle, the 
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singular behavior observed could disappear in the ther- 
modynamic limit, however, we find that this is not the 
case. 

By definition, the threshold temperature is a function 
of parameters (ft, rj) and the system size N, T t h(h, rj : N). 
In this subsection, we show that the threshold tempera- 
ture for large N is independent of the system size N, 
namely, T t h(h,r) : N) = T t h(h,rj) for large N. It is 
numerically confirmed that N = 80 is large enough for 
achieving the thermodynamic limit for this system. 

The proof of T th (h, rj : N) — T th (h, rj) is based on the 
following two facts. First, the geometric measure of the 
ground state of the XY model scales polynomially with 
N for large N, G(po) oc C(ft, rj)N where a function ((h, rj) 
is determined independently of N. Second, in large sys- 
tems, for any fi G [0, 1] independent of N, there exists an 
unique temperature T{ji) such that po(h, rj : T(fi) : N) — 
2~ftN _ rp Q that, we prove that the scaling of the pop- 
ulation of the ground state Po(h, rj : T : N) is exponential 
with the system size N, p (h,ij : T : N) cx 2~^ h ^ N 
where log 2 (l + e'^+^Z ^ bT '>) < fj,(h,rj : T) < log 2 (l + 
e -8(K n )/(k B T)^ where d(h,r/) = \h- 1| for ft > 1 -rj 2 and 

d(h,rj) = T)\/l — if- — b?/yj\ — rj 2 for h < 1 — rj 2 . Since 
the exponent fJb{h, rj : T) is monotonically increasing with 
temperature, the statement holds. 

From the first fact, the definition of the threshold tem- 
perature, po(h,rj : T t h '■ N) = 2~ G ( po \ is rewritten as 



Po 



(h,r,:T th (h,r l :N):N) = 2-^ N (16) 



for large N. On the other hand, due to the second fact, 
there exists an unique temperature T(£(/i, 77)) such that 



Po 



(M : T(C(M)) : AO = 2" C(/M,) 



N 



(17) 



if N is large enough. Due to these two equations and 
the fact that the population of the ground state is a 
single-valued function with respect to temperature, we 
can conclude that T t h(h, rj : N) — T(£(h, rj)) in large sys- 
tems, namely the threshold temperature is independent 
of N for large N. 



V. THERMAL ROBUSTNESS OF 
ENTANGLEMENT OF THE GROUND STATE 

We now look in more detail at the value of the thresh- 
old temperature, beyond the singular behavior of the 
threshold temperature. In the study of the threshold 
temperature related to the phase diagram at zero tem- 
perature, it is enough to analyze its singularity since the 
phase diagram is also determined by the singular behav- 
ior of the correlation functions, in which quantitative dif- 
ferences are negligible. 

On the other hand, we can study the threshold temper- 
ature as an indicator of the thermal robustness of entan- 
glement of the ground state by investigating its value. 
It may be natural to think that the stronger (weaker) 
the entanglement of the ground state, the higher (lower) 



the temperature at which the entanglement can survive 
in a thermal state. However, this is not always correct. 
In this section, we describe the characteristic behaviors 
of multipartite entanglement caused by the existence of 
thermal excitations in two regions, one is on the H-line, 
and the other is for ft> 1. In those regions, we have find 
this natural intuition is broken. 



A. Threshold temperature on the H-line 

In Fig. [3] the geometric measure of the ground state 
takes the maximum on the H-line. This fact results in 
the local maximum of the threshold temperature on the 



line, as we explained in Subsection IV A But the peak of 
the threshold temperature on the line is less emphasized 
compared to that of the geometric measure. As long as 
the population of the ground state does not change so 
much with parameters, the maximum of the geometric 
measure is inherited by the threshold temperature, and 
the peak of the maximum of the threshold temperature 
should be more pronounced than that of the geometric 
measure. However, on the H-linc, the population of the 
ground state takes local minimum since the first excited 
energy gap, between the ground energy and the first ex- 
cited energy, approaches zero and the system reduces to 
a gapless system in the thermodynamic limit. This local 
minimum of the population of the ground state leads a 
smoother peak of the threshold temperature on the line. 

Therefore, there is a balance between the enhancement 
of the peak due to the local maximum of entanglement 
of the ground state and the smoothing of the peak due to 
the local minimum of the population of the ground state. 
As a net result, the peak of the threshold temperature 
is not enhanced much. Thus, multipartite entanglement 
at zero temperature is high on the H-line, but is fragile 
against thermal excitations. 



B. Threshold temperature in the region of h ^> 1 

We now look at the threshold temperature in the XY 
model for h ^> 1. On the line of h > 1 and 77 = the 
ground state is separable and, in the region of ft, > 1 
and rj ^ 0, the geometric measure of the entanglement of 
the ground state has small but nonzero value due to the 
effect of nonzero rj. The threshold temperature is given 
as in Fig. [6] We present only for the case of 77 € [0,1]. 

We can see that the threshold temperature increases 
monotonically with the external magnetic field h in the 
region of ft ^ 1 for rj ^ 0. Since a thermal state below 
the threshold temperature is guaranteed to be entangled, 
a thermal state of the XY model for ft ^ 1 is entangled 
even at high temperature, although the geometric mea- 
sure of the ground state is monotonically decreasing in 
that region. It may be counterintuitive that such weak 
entanglement of the ground state can survive in a thermal 
state at high temperature. The high threshold temper- 
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1 1_ ( to quantum properties at the limit of h 3> 1, it is interest- 

ing to investigate how multipartite entanglement at finite 
temperature contributes to quantum properties of the 
„-*■"' system in the trade-off of thermal robustness of entan- 

.^-'** glcmcnt and the amount of entanglement of the ground 

," m " state in the intermediate region h > 1. 



fl". h VI. CONCLUSION 

10 20 30 40 50 



FIG. 6: (Color online) The threshold temperature in the XY 
model calculated in the system with N = 80. We set Boltz- 
mann constant kn as one. The (red) solid line with small 
squares, the (orange) dashed line, the (green) dotted dashed 
line and the (blue) dotted line show the threshold tempera- 
ture for 77 = 0,0.1,0.4,0.9, respectively. The parameter h is 
between and 50. We can see that the threshold temperature 
is monotonically increasing with h in the region of h 3> 1. We 
note that the threshold temperature for the case of r\ — 
and ft > 1 is zero since the ground state is separable in that 
case. We can see that the effect of nonzero rj to the threshold 
temperature is high although entanglement itself is small (see 
Fig. 

ature originates from the fact that the first energy gap 
between the ground energy and the first excited energy 
increases more rapidly than the rate of decrease of the ge- 
ometric measure. Thus weak entanglement of the ground 
state does not necessarily imply that it is easily destroyed 
by thermal excitations. 

We note that the high threshold temperature does not 
mean thermal effects can boost the amount of the entan- 
glement. The ground state in the region of h ^> 1 for 
i] is very close to a product state and the amount of 
entanglement is infinitcsimally small but entanglement is 
thermally robust in a sense that it can survive until high 
temperature. In relatively artificial systems considered 
in quantum information, it has been found [28) that a 
multipartite entangled state which is close to a product 
state can be useful as a universal resource for the mea- 
surement based quantum computation |29j . In this case, 
the existence of infinitesimal entanglement provides a sig- 
nificant difference from zero-entanglement cases. Also, 
a related approach on entanglement to that presented 
here provides an entanglement witness very closely linked 
to the threshold temperature which manifest themselves 
as physical observable in some spin systems |30j . How- 
ever, so far, we do not know whether such an infinitesi- 
mal amount but thermally robust entanglement can con- 
tribute to present any quantum property of natural phys- 
ical systems at finite temperature. It is desirable to inves- 
tigate physical properties at finite temperature related to 
entanglement of the ground state in the region of h 3> 1 . 
On the other hand, even if it would turn out that the in- 
finitesimal amount of entanglement does not contribute 



In this paper, we have studied thermal robustness of 
multipartite entanglement of the ground state of the 
one-dimensional spin 1/2 XY model with a transverse 
magnetic field by using the threshold temperature in- 
troduced by |17j . The threshold temperature is a tem- 
perature where the thermal state of a multipartite sys- 
tem is certainly entangled and its lower bound of the 
temperature can be obtained in terms of the geometric 
measure. We have calculated the geometric measure and 
threshold temperature in the spin 1/2 one dimensional 
XY model. We have analyzed the correspondence be- 
tween the threshold temperature and the phase diagram 
at zero temperature by looking at the singularity of the 
threshold temperature. We have also analyzed thermal 
robustness of multipartite entanglement of the ground 
state. 

It has been shown that the C-line, on which two point 
correlation functions are constant and entanglement of 
the ground state takes local minimum with respect to 
parameters in the Hamiltonian, is clearly reflected by 
the threshold temperature as well as the V-line (the sec- 
ond order QPT line) and the H-line (the first order QPT 
line). Although the geometric measure of the ground 
state itself does not reflect the C-linc clearly, the C-line 
is enhanced by taking into account thermal excitations. 
We also present a region in the phase diagram where 
multipartite entanglement at zero temperature is high 
but is thermally fragile, and another region where mul- 
tipartite entanglement at zero temperature is less but is 
thermally robust. These entanglement properties can- 
not be deduced only from analyzing the ground state. 
They provide novel information of the quantum effects 
corresponding to QPT for natural physical systems at 
the non-zero temperature. 
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APPENDIX A: GAPPED THRESHOLD 
TEMPERATURE 



APPENDIX B: EIGENSTATES AND 
EIGENENERGY OF THE XY MODEL 



We show that 

(l_ e -A/(*3T) )po(T)> 



1 



d [R(p ) + l] 
p t h(T) is entangled. (Al) 

Let pi := Pi/di be a mixed state for the i-th excitation. 
A thermal state p t h{T) can be rewritten as 

Pth(T) = (po(T)d +p 1 (T)d 1 ) 

1 



x [poi + ( 



Po(T)d + pi(T)di 



l)p r 



where p,(T) = exp[-Ei/(k B T)]/Z(T) and 
p (T)d po + Pi{T)dxpi 



Pox 



Prest 



Po(T)d 



■Pi(T)di 
1 



1 - ( Po (T)d 



PidrPi 



(A2) 

(A3) 
(A4) 



By definition of the robustness of entanglement, 

1 < R(poi) => Pth(T) is entangled. 



Po(T)d +p 1 (T)d 1 

(A5) 

We relate the robustness of entanglement of the state poi, 
R(pox), to that of the ground state, R(po). The state poi 
can be divided into the ground state po an d the another 
state as 

e -A/(k B T) )do 

-Po 



Poi 



(1 

d + e^ A /( fe « T )di 
(d + rfi)e- A/(fcBT) d oPo + d lPl 



(A6) 



d + e- A /( fc « T )di d + di 

Since the robustness of entanglement satisfies 

R(pa+(l-p)a') + l >p[R(a) + l] (A7) 

for p > 0, p € [0, 1] and any states <r, a', the robustness 
R(pox) and R(po) can be related as 



#(poi) + 1 > 



(1 



-A/(k B T) 



[R( Po ) + l]. (AS) 



d + e -A/(fc B r) dl 

By replacing the right-hand side of the inequality in 
Eq. ( A5 1 , we can obtain the inequality that 

1 



( l_ e -A/( feB T) )po(T)> 



(A9) 



d [R( Po ) + iy 

The reason why the gapped threshold temperature is 
always smaller than the threshold temperature is that 



the bound of the robustness of entanglement ( A7 ) is not 



tight. In general, the threshold temperature is expected 
to be improved by taking into account of the first energy 
gap as we can see from the proof in Ref. [T7] . 

We also note that it is possible to bound po simply 
in terms of the gap A, as was done in [17], however, in 
general this is not useful since it gives very low threshold 
temperatures, going to zero in the thermodynamic limit. 



We present a brief explanation of the XY model. The 
Hamiltonian is given by 



H 



XY 



N 

£ 

j=x 



^~ a 3 ff J+l 



1 



(Bl) 

Although we consider -q S [—1,1], it is enough to study 
the case of r\ € [0, 1]. This is because the Hamiltonian 

in the region of 77 6 [—1,0), say H X y , can be obtained 
by multiplying a local unitary by the Hamiltonian in the 
region of rj 6 (0, 1], say H XY \ namely 



H 



[-1,0) _ 



XY 



(B2) 



where Ui — (I + iaf) / '\p2. Hence we can easily construct 
the eigenstates and the eigenenergy of the Hamiltonian 
from those of the Hamiltonian H XY ' ■ Similarly, 



H 



-1,0) 



XY 

we can extend the region of parameters to arbitrary h 
by multiplying by another local unitary. In the following 
explanation of the XY model, we consider only the case 
of 77 G [0,1]. 

This Hamiltonian is exactly diagonalized by the 
Jordan- Wigner transformation Ck '■— ®jZ\(—af)a^ 
where a k := (a k ± ia~^)/\/2, the Fourier transforma- 
tion f k := ^=Ejli ex P[-(0fe j + j)i]cj and the Bo- 



goliubov transformation T k (h,rj) 



sin <p k (h, r])f J, where 9 
{4>^{h, r])} are given by 



5+ 



cos (^{h,r])f k - 

e*=ik a- - »i and 



N 



N 



rj sin 9 k cos 2<jy^ 
r\ sin 9 k sin 2(jr k 



-(h- 



0. 



cos 9 k ) sin 
s6%) cos 2^ > 



i2</>± 



cos 



(B3) 
(B4) 



Since the Hamiltonian Hxy commutes with parity op- 
erator — 2a^cr~), we can divide the whole Hilbcrt 

space TixY into a parity plus (minus) subspace, 7~i XY , as 



'Hxy = H 



XY 



n 



XY' 



(B5) 



The states in a parity plus (minus) subspace have even 
(odd) number of up spins. The Hamiltonians H XY pro- 
jected onto each subspace Ti. XY are given by 



N 1 
Hxy =£4(M)(lf f (M)lt(M - x), 



(B6) 



k=l 



where e ± (k) := 2^(h + cos 6^) 2 + r) 2 sin 2 0±(> 0). 

Since the operators V k ' (k = 1, • • • ,N) satisfy the anti- 
commutation relation, they can be regarded as genera- 
tion operators of pseudo fermions with an energy e (k). 
The lowest energy for each Hamiltonian, H XY or H XY , 
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is given by 



E+(h,rj) = -^^/(/i + cosfl+^ + T^sin 2 



'k ' 



(B7) 



k=l 



N 



e (h, 77) = -E(h) -J2v( h+ cos 



fc=i 



9,7) 2 + ?7 2 sin 2 ^, 
(B8) 



respectively, where = for /i < 1 and = 2(h — 
1) for /i > 1. The eigenstates corresponding to the lowest 
energy in each subspaces are given by 



|e+) = n(cofl^++sin^+/+V^)| I •■■ I), (B9) 

k=l 



fc=i 



(BIO) 



In the region of h > 1, the ground energy E {h,rj) is 
equal to E^ih,^) and the ground state is |ej). In the 
region of h < 1, two energies E^(h,rj) coincide in the 
thermodynamic limit (iV — > 00) and the ground states 
degenerate. Since a thermal state at zero temperature 
reduces to a projector onto the subspace spanned by the 
degenerated ground state, the ground state in the ther- 
modynamic limit in the region of h < 1 is a mixed state. 

The eigenstates |eg) in Eq. (B9l and Eq. (BIO I are 
represented by global operators {f^}k- To investigate 
entanglement of the states, the states should be repre- 
sented in a local basis. Thus we have to transform global 
operators back to local operators defined at each site. 



APPENDIX C: INVESTIGATION OF (^(M) ON 
THE V-LINE AND ON THE H-LINE 

We show that (4(1 - Sh,rf)\e^{l + 8h,rf)) = for 
any small Sh > in the thermodynamic limit, which had 
been numerically obtained in Ref. [27] For that, we first 



investigate the varia bles {(f> k [h, rj)}k defined by Eq. (B3) 



and the inequality (B4|, which are rewritten as 



V 



sm6f 



tan 20± (h, 77) = - — . 

h + cos 0jjr 

i20±(M) >o. 



sin 1 



(CI) 
(C2) 



Since 62 



>k ~ (2fe ^' K = % the variables {<ft(h,p)} k 
reduces to {4>e(h, i])}oe(o,Tr)i in the thermodynamic limit, 
defined by 



tan 2<pe(h, rj) = — 
sin 20 e (M) > 0. 



r) sin 9 
h + cos 9 ' 



(C3) 
(C4) 




IT/2 

3-/8 






■* -t/4 






jt/8 








-/4 -12 3-/4 - 







~1tA TT2t~ ~ 
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-0.5 
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np. 


-1 


3-/8 


-1.25 


* .t/4 
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<r/B 


-1.75 




-2 





FIG. 7: (Color online) The upper figure is the numerical 
calculations of 4>„ /2 (h, 0.5) for N = 10, 20, 50, 100, 1000. We 
can see the non-analyticity at h = 1. The middle figure and 
the lower figure are the right-hand side of Eq. ( C3 1 in the 
case of h = 0.9 (middle) and h ■ 
6. We fixed 77 = 0.5. 
respectively. 



1.1 (lower) with respect to 
The inset is <f>e (h, 77) with respect to 9 



We show that, for any small Sh > 0, 
lim 4>e(l — Sh, rj) = 0, 

7T 

lim 4> e {l + Sh,rf) = - 



(C5) 
(C6) 



In the case of h = 1 — Sh, there exists 9q such that the 
right-hand side of Eq. ( C3 1 diverges. Then tan20e(l — 
Sh, 77) > for arccos(— h) < 9 < n (see Fig. [7]). Therefore 
lime_, ff 0e(l — Sh,r)) = 0. On the other hand, in the 
case of h = 1 + Sh, the right-hand side of Eq. ( C3 1 is 
always finite, which implies tan20e(l — Sh,r/) < for any 
9 G (0,tt) (see Fig. M. Then lim^ <f> g (1 + Sh, 77) = f . 
Due to this fact, we can easily show that (ej(l — 



*M)|eo(l + *M)> 
thermodynamic limit 



= for any small Sh > in the 
The inner product can be rewrit- 
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ten as 



(4{1-Sh,rf)\4(l + Sh,v)) 



M 



= Y[cos[^(l + 5h)-4>t(l-Sh)], (C7) 

k=l 

where AT is N/2 for the plus case and N/2 — 1 for 
the minus case. Since Eq. (C5l and Eq. (C6l imply 
(j>Af{l + Sh,r]) — 4>^f(l — Sh, rj) = | in the thermody- 
namic limit, the inner product is zero. As long as no 
accidental cancelation of terms occurs in the calculation 



of the geometric measure, the geometric measure inher- 
its this abrupt change of the ground state and has non- 
analyticity on the V-line as explained in Subsection |III B| 



In a similar 



way, it is proven that 
(e^(/i, 5r))\e^(h, — Srf)) — for any small 8r\ > in 
the thermodynamic limit. Then the ground state also 
has an dramatic change around the H-line. 

We note that, since the variables {4>g(h, r/)} completely 
determine how the ground state depends on the parame- 
ters (h,rf), the non-analyticity of these variables will in- 
herit to the ground state. 
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